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On the way to meet the experimental observation of persistent current in a 
mesoscopic cyhnder: A mean field study 
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The behavior of persistent current in a mesoscopic cylinder threaded by an Aharonov-Bohm flux 
<j) is carefully investigated within a Hartree-Fock mean field approach. We examine the combined 
effect of second-neighbor hopping integral and Hubbard correlation on the enhancement of persistent 
current in presence of disorder. A significant change in current amplitude is observed compared to 
the traditional nearest-neighbor hopping model and the current amplitude becomes quite comparable 
to experimental realizations. Our analysis is found to exhibit several interesting results which have 
so far remained unaddressed. 
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I. INTRODUCTION 



Over the last many years appearance of persistent cur- 
rent in metallic single-channel rings and multi-channel 
cylinders has drawn much attention in theoretical as well 
as in experimental research. In mesoscopic regime where 
dimensions of a system are comparable to mean free path 
of an electron the phase coherence of electronic states is of 
fundamental importance and the existence of dissipation- 
less current in a mesoscopic conducting ring threaded by 
an Aharonov-Bohm (AB) flux is a direct consequence of 
quantum phase coherence. In this new quantum regime, 
two important aspects appear at low temperatures and 
they are as follows. 

• The phase coherence length L^ i.e., the length scale for 
which an electron maintains its phase memory, increases 
significantly with the lowering of temperature and be- 
comes comparable to the system size L. 

• The energy levels of these small finite size systems are 
discrete. 

These two are the most essential criteria for the existence 
of persistent current in a small metallic ring/cylinder due 
to the application of an external magnetic flux </>. In 
the pioneering work of Biittiker, Imry and Landauer [l|, 
the appearance of persistent current in metallic rings has 
been explored. Later, many excellent experiments |2t2| 
have been carried out in several ring and cylindrical ge- 
ometries to reveal the actual mechanisms of persistent 
current. Though much efforts have been paid to study 
persistent current both theoretically ]&h2^ as well as ex- 
perimentally 0-01, yet several drawbacks still exist be- 
tween the theory and experiment, and the full knowledge 
about it in this scale is not well established even today. 

The results of the single loop experiments are signif- 
icantly different from those for the ensemble of isolated 
loops. Persistent currents with expected (/)o periodicity 
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have been observed in isolated single Au rings Q and 
in a GaAs-AlGaAs ring [j|. Levy et al. [2| found oscil- 
lations with period 0o/2 rather than (f>Q in an ensemble 
of 10^ independent Cu rings. Similar (f>o/2 oscillations 
were also reported for an ensemble of disconnected 10^ 
Ag rings Q as well as for an array of 10^ isolated GaAs- 
AlGaAs rings [7|. In a recent experiment, Jariwala et 
al. [5| obtained both (po and (/f)o/2 periodic persistent cur- 
rents in an array of thirty diffusive mesoscopic Au rings. 
Except for the case of the nearly ballistic GaAs-AlGaAs 
ring [j|, all the measured currents are in general one or 
two orders of magnitude larger than those expected from 
the theory. 

Free electron theory predicts that at T = 0, an ordered 
one-dimensional metallic ring threaded by magnetic flux 
(j) supports persistent current with maximum amplitude 
/o = evp / L, where vp is the Fermi velocity and L is the 
circumference of the ring. Metals are intrinsically dis- 
ordered which tends to decrease the persistent current, 
and the calculations show that the disorder-averaged cur- 
rent (/)crucially depends on the choice of the ensem- 
ble [I, [il, [il. The magnitude of the current (/2)i/2 jg 
however insensitive to the averaging issues, and is of the 
order of IqI/L, I being the elastic mean free path of the 
electrons. This expression remains valid even if one takes 
into account the finite width of the ring by adding con- 
tributions from the transverse channels, since disorder 
leads to a compensation between the channels [3, [l3| . 
However, the measurements on an ensemble of 10^ Cu 
rings [2| reported a diamagnetic persistent current of av- 
erage amplitude 3 x 10""^ evp/L with half a flux-quantum 
periodicity. Such (j)o/2 oscillations with diamagnetic re- 
sponse were also found in other persistent current exper- 
iments consisting of ensemble of isolated rings [g, 01 • 

Measurements on single isolated mesoscopic rings on 
the other hand detected (/)o-periodic persistent currents 
with amplitudes of the order of Iq ~ evp/L, (closed to 
the value for an ordered ring). Theory and experiment [J] 
seem to agree only when disorder is weak. In another re- 
cent nice experiment Bluhm et al. (30| have measured 
the magnetic response of 33 individual cold mesoscopic 



gold rings, one ring at a time, using a scanning SQUID 
technique. They have measured h/e component and pre- 
dicted that the measured current amphtude agrees quite 
weh with theory [g| in a single ballistic ring [J] and an en- 
semble of 16 nearly ballistic rings [31[. However, the am- 
plitudes of the currents in single-isolated-diffusive gold 
rings [3| were two orders of magnitude larger than the 
theoretical estimates. This discrepancy initiated intense 
theoretical activity, and it is generally believed that the 
electron-electron correlation plays an important role in 
the disordered diffusive rings |19l42l| . An explanation 
based on the perturbative calculation in presence of in- 
teraction and disorder has been proposed and it seems 
to give a quantitative estimate closer to the experimen- 
tal results, but still it is less than the measured currents 
by an order of magnitude, and the interaction param- 
eter used in the theory is not well understood physi- 
cally. Most of these theoretical results have been ob- 
tained based on a tight-binding framework within the 
nearest-neighbor hopping (NNH) approximation. This is 
an important approximation and it has been shown that 
within the NNH model electronic correlation provides a 
small enhancement of current amplitude in disordered 
materials i.e., a weak delocalizing effect is observed in 
presence of electron-electron (e-e) interaction. As an at- 
tempt in the present work we modify the traditional NNH 
model by incorporating the effects of higher order hop- 
ping integrals, at least second- neighbor hopping (SNH), 
in addition to the NNH integral. It is also quite physical 
since electrons have some finite probabilities to hop from 
one site to other sites apart from nearest-neighbor with 
reduced strengths. We will show that the inclusion of 
higher order hopping integrals gives significant enhance- 
ment of current amplitude and it reaches quite closer to 
the current amplitude of ordered systems. 

The other important controversy comes for the deter- 
mination of the sign of low-field currents and still it is 
an unresolved issue between theoretical and experimental 
results. In an experiment on persistent current Levy et 
al. [4I have shown diamagnetic nature for the measured 
currents at low-field limit. While, in other experiment 
Chandrasekhar et al. 3 have obtained paramagnetic re- 
sponse near zero field limit. Jariwala et al. [5| have pre- 
dicted diamagnetic persistent current in their experiment 
and similar diamagnetic response in the vicinity of zero 
field limit were also supported in an experiment done by 
Deblock [y] et al. on Ag rings. Yu and Fowler [22l have 
shown both diamagnetic and paramagnetic responses in 
mesoscopic Hubbard rings. Though in a theoretical work 
Cheung et al. [^ have predicted that the direction of cur- 
rent is random depending on the total number of elec- 
trons in the system and the specific realization of the 
random potentials. Hence, prediction of the sign of low- 
field currents is still an open challenge and further studies 
on persistent current in mesoscopic systems are needed 
to remove the existing controversies. 

In the present paper we address the behavior of persis- 
tent current in an interacting mesoscopic ring with finite 



width threaded by an Aharonov-Bohm flux <f). A simple 
tight-binding Hamiltonian is used to illustrate the sys- 
tem and all the calculations are performed within a mean 
field approach. Using a generalized Hartree-Fock (HF) 
approximation |32l . |33| , we compute numerically persis- 
tent current (/) as functions of AB fiux (/), total number of 
electrons N^.^ e-e interaction strength U , second- neighbor 
hopping integral, disorder strength W and system size N . 
The main motivation of the present work is to illustrate 
the effects of higher order hopping integrals on the en- 
hancement of persistent current in disordered mesoscopic 
cylinders. Our results can be utilized to study magnetic 
response in any interacting mesoscopic system. 

In what follows, we present the results. In Section II, 
we describe the geometric model and generalized Hartree- 
Fock theory to study magnetic response in the model 
quantum system. Section HI contains the numerical re- 
sults. Finally, in Section IV we draw our conclusions. 



II. MODEL AND SYNOPSIS OF THE 
THEORETICAL FORMULATION 



Let us start by referring to Fig.[TJ where a small metal- 
lic cylinder is threaded by a magnetic fiux cj). The filled 
black circles correspond to the positions of the atomic 
sites in the cylinder. To predict the size of a cylinder we 




FIG. 1: (Color online). Schematic view of a ID mesoscopic 
cylinder penetrated by a magnetic flux (p. The red dashed 
line corresponds to the second-neighbor hopping integral and 
the filled black circles represent the positions of the atomic 
sites. A persistent current / is established in the cylinder. 



use two parameters N and M, where the 1st one (N) rep- 
resents total number of atomic sites in each circular ring 
and the other one (M) gives total number of identical 
circular rings. For the description of our model quantum 
system we use a tight-binding (T-B) framework and in 
order to incorporate the effect of higher order hopping 
integrals to the Hamiltonian here we consider second- 
neighbor hopping (SNH) (shown by the red dashed line 
in Fig. [1} in addition to the nearest-neighbor hopping 
(NNH) of electrons. Considering both NNH and SNH 
integrals the T-B Hamiltonian for the cylindrical system 



in Wannicr basis looks in the form, 



written in the decoupled form as, 



He = Xl^''J''^4j,'TCj,i,<T + Xl^' ['='^'4j,<t'=*J+i^'^ 



h.c] + '^td le^'^'^cl ,,■ ,,Q+i,j+i,<j + h.c. 



(1) 



where, (?,j) represent the co-ordinate of a lattice site. 
The index i runs from 1 to M, while the integer j goes 
from 1 to N. eij-,o- is the on-site energy of an electron 
at the site {i,j) of spin a (t,i)- ti and td are the NNH 
and SNH integrals, respectively. Due to the presence of 
magnetic flux (f> (measured in unit of the elementary flux 
quantum 0o = ch/e), a phase factor 9i = 2Tr(j)/N appears 
in the Hamiltonian when an electron hops longitudinally 
from one site to its neighboring site, and accordingly, a 
negative sign comes when the electron hops in the re- 
verse direction. 6d is the associated phase factor for the 
diagonal motion of an electron between two neighboring 
concentric rings. No phase factor appears when an elec- 
tron moves along the vertical direction which is set by 
proper choice of the gauge for the vector potential A as- 
sociated with the magnetic field B, and this choice makes 
the phase factors {6i, 6d) identical to each other for the 
longitudinal and diagonal motions. Since the magnetic 
field corresponding to the AB flux does not penetrate 
anywhere of the surface of the cylinder, we ignore Zeeman 
term in the above tight-binding Hamiltonian (Eq. [ij. 
c]^ g. and Ci.j^a are the creation and annihilation oper- 
ators, respectively, of an electron at the site (i,j) with 
spin a. U \s the on-site Hubbard interaction term. 



A. Decoupling of the interacting Hamiltonian 



Hr 



i,j i-'j 

Ytd[e:^'c\j,.c,+i,j+i,^ + h.c. 

i-j 

^t + ^i-E^("^J^t>KM) (4) 



^,3 



where, H^ and H^ correspond to the effective tight- 
binding Hamiltonians for the up and down spin electrons, 
respectively. The last term is a constant term which pro- 
vides an energy shift in the total energy. 



B. Self consistent procedure 

With these decoupled Hamiltonians {H^ and H^) of 
up and down spin electrons, now we start our self consis- 
tent procedure considering initial guess values of (j^ij-.t) 
and {ni,j^^). For these initial set of values of {nij,-^) and 
('^ij,^), we numerically diagonalize the up and down spin 
Hamiltonians. Then we calculate a new set of values of 



'■'1,3 



f) and {ni,j^i). These steps are repeated until a self 



consistent solution is achieved. 



To get the energy eigenvalues of the interacting model 
quantum system described by the above tight-binding 
Hamiltonian given in Eq. [U first we decouple the inter- 
acting Hamiltonian using generalized Hartree-Fock ap- 
proach, the so-called mean field approximation. In this 
procedure, the full Hamiltonian is completely decoupled 
into two parts. One is associated with the up-spin elec- 
trons, while the other is related to the down-spin elec- 
trons with their modified site energies. For up and down 
spin Hamiltonians, the modified site energies are ex- 
pressed in the form. 



^i,j:t 



e.,j,t + U{nij^i} 



(2) 



C. Calculation of ground state energy 

After achieving the self consistent solution, the ground 
state energy Eg for a particular filling at absolute zero 
temperature (T = OK) can be determined by taking the 
sum of individual states up to Fermi energy (Ep) for both 
up and down spins. Thus, we can write the final form of 
ground state energy as, 

E0=Y1 ^P^t + Y ^P-i - Y ^("*.^^t> irl^.3,i) (5) 



4j4 = "^^^J.^ + U{ni,j,t) 



(3) 



where. 



''J.o- 



'^i,j,a-'^'',3,cr 



is the number operator. With 



these site energies, the full Hamiltonian (Eq. [T]) can be 



where, the index p runs for the states up to the Fermi 
level. Ep^^ (Ep^i) is the single particle energy eigenvalue 
for p-th eigenstate obtained by diagonalizing the Hamil- 
tonian H^ (Hi). 



D. Calculation of persistent current 

At absolute zero temperature, total persistent current 
of the system is obtained from the expression, 



/(0) = -c 



dE^icp) 



(6) 



where, Eq{(J)) is the ground state energy for a particular 
filling. 

In the present work we perforin all the essential fea- 
tures of persistent current at absolute zero temperature 
and use the units where c = h = e = 1. Throughout our 
numerical calculations we set the nearest-neighbor hop- 
ping strength ti = —1 and fix M = 2 i.e., cylinders with 
two identical rings. Energy scale is measured in unit of 
ti. 
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III. NUMERICAL RESULTS AND DISCUSSION 

Following the above theoretical prescription now we 
start to analyze our numerical results. We describe the 
results in three different parts. In the first part, we con- 
sider perfect cylinders with only nearest-neighbor hop- 
ping integral. In the second part, disordered cylinders 
described with only NNH integral are considered. Fi- 
nally, in the third part we discuss the effect of second- 
neighbor hopping (SNH) integral on the enhancement of 
persistent current in disordered cylinders. 



A. Perfect cylinders with NNH integral 

For perfect cylinders we choose £i,j.-\- = ^i,j,i = for all 
(i, j). Since here we consider the cylinders described with 
NNH integral only, the second-neighbor hopping strength 
td is fixed to zero. 



FIG. 2: (Color online). Ground state energy levels as a func- 
tion of flux (j) for some perfect cylinders with A^ = 5 and 
M — 2. The red, green and blue curves correspond to U — 0, 
0.5 and 1, respectively, (a) Quarter-filled case and (b) Half- 
flUed case. 
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Energy-flux characteristics 




As illustrative examples, in Fig. [2] we show the varia- 
tion of ground state energy levels as a function of mag- 
netic fiux (j) for some typical mesoscopic cylinders where 
N is fixed at 5 (odd N). In (a) the results are given 
for the quarterly-filled [N^ = 5) cylinders, while in (b) 
the curves correspond to the results for the half-filled 
[Nf, — 10) cylinders. The red, green and blue lines rep- 
resent the ground state energy levels for t/ = 0, 0.5 and 
1, respectively. It is observed that the ground state en- 
ergy shows oscillatory behavior as a function of cf) and the 
energy increases as the electronic correlation strength U 
gets increased. Most significantly we see that the ground 
state energy levels provide two different types of peri- 
odicities depending on the electron filling. At quarter- 
filling, ground state energy level gives 4>o (= 1, since 
c = e = /i=lin our chosen unit system) flux-quantum 



FIG. 3: (Color online). Ground state energy levels as a func- 
tion of flux <j) for some perfect cylinders considering N — 8 
and M = 2. The red, green and blue curves correspond to 
U = 0, 0.5 and 1, respectively, (a) Quarter-flUed case and (b) 
Half-filled case. 



periodicity. On the other hand, at half-filling it shows 
(/>o/2 flux-quantum periodicity. The situation becomes 
quite different when the total number of atomic sites N 
in individual rings is even. For our illustrative purposes 
in Fig. [3] we plot the lowest energy levels as a function 
of (f) for some typical mesoscopic cylinders considering 
A^ = 8 (even N). The curves of different colors corre- 
spond to the identical meaning as in Fig. [21 From the 
spectra given in Figs. [2Ia) (quarter-flllcd case) and (b) 



(half-filled case) it is clearly observed that the ground 
state energy levels vary periodically with AB flux (jj ex- 
hibiting only 00 flux-quantum periodicity. Thus it can 
be emphasized that the appearance of half flux-quantum 
periodicity strongly depends on the electron filling as well 
as on the oddness and evenness of the total number of 
atomic sites N in individual rings. Only for the half-filled 
cylinders with odd N, the lowest energy level gets 0o/2 
periodicity with flux (j). Now it is important to note that 
this half flux-quantum periodicity does not depend on the 
width (M) of the cylinder and also it is independent of 
the Hubbard correlation strength U. Hence, depending 
on the system size and filling of electrons variable peri- 
odicities are observed in the variation of lowest energy 
level. It may provide an important signature in studying 
magnetic response in nano-scale loop geometries. 



2. Current-flux characteristics 

In Fig. m we display the current-flux characteristics 
for some impurity free mesoscopic cylinders considering 
M ~ 2. In (a) the results are given for the half- filled case 
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FIG. 4: (Color online). Persistent current as a function of flux 
(j> for some ordered mesoscopic cylinders considering M = 2. 
(a) Half-filled case with A*' — 15. The red, green and blue 
curves correspond to U — 0, 1.5 and 2, respectively, (b) 
Quarter-filled case with A*' — 20. The red, green and blue 
curves correspond to [/ = 0, 2 and 3, respectively. 

where we set N ~ 15. The red line corresponds to the 
current for the non-interacting {U = 0) case, while the 
green and blue lines represent the currents when U = 1.5 
and 2, respectively. From the curves we notice that the 
current amplitude gradually decreases with the increase 
of electronic correlation strength U. The reason is that 
at half-filling each site is occupied by at least one elec- 
tron of up spin or down spin, and the placing of a second 
electron of opposite spin needs more energy due to the 




repulsive effect of U. Thus conduction becomes difficult 
as it requires more energy when an electron hops from 
its own site and situates at the neighboring site. Now 
both for the non-interacting and interacting cases, cur- 
rent shows half flux-quantum periodicity as a function a 
(j) obeying the energy-flux characteristics since here we 
choose odd N {N = 15). The behavior of the persis- 
tent currents for even N is shown in (b) where we set 
A^ = 20. The currents are drawn for the quarter-filled 
case i.e., N^ = 20, where the red, green and blue curves 
correspond to [/ = 0, 2 and 3, respectively. The reduc- 
tion of current amplitude with the increase of Hubbard 
interaction strength is also observed for this quarter-filled 
case, similar to the case of half-filled as described earlier. 
But the point is that at quarter-filling, the reduction of 
current amplitude is much smaller compared to the half- 
filled situation. This is quite obvious in the sense that 
at less than half-filling 'empty' lattice sites are available 
where electrons can hop easily without any cost of ex- 
tra energy and the conduction becomes much easier than 
the half-filled situation. In this quarter-filled case, per- 
sistent currents provide only 0o flux-quantum periodicity 
following the E-cj) diagram. From these current-flux char- 
acteristics it can be concluded that for 'ordered ' cylinders 
current amplitude always decreases with the enhancement 
in Hubbard correlation strength U . 



B. Disordered cylinders with NNH integral 

In order to describe the effect of impurities on elec- 
tron transport now we focus our attention on the results 
of some typical disordered cylinders described with NNH 
integral. Here we consider the diagonal disordered cylin- 
ders i.e., impurities are introduced only at the site ener- 
gies without disturbing the hopping integrals. The site 
energies in each concentric ring are chosen from a corre- 
lated distribution function which looks in the form. 



ej,t = ej.i = W^cos(jA7r) 



(7) 



where, W is the impurity strength. A is an irrational 
number and we choose A = {l + ^/5)/2, for the sake of our 
illustration. Setting A = 0, we get back the pure system 
with uniform site energy W. Now, instead of considering 
site energies from a correlated distribution function, as 
mentioned above in Eq. [71 we can also take them ran- 
domly from a "Box" distribution function of width W. 
But in the later case we have to take the average over 
a large number of disordered configurations (from the 
stand point of statistical average) and since it is really a 
difficult task in the aspect of numerical computation we 
select the other option. Not only that in the averaging 
process several mesoscopic phenomena may disappear. 
Therefore, the averaging process is an important issue in 
low-dimensional systems. 

In presence of disorder, energy levels get modified sig- 
nificantly. For our illustrative purposes in Fig. [5] we plot 
ground state energy levels as a function of magnetic fiux 



(j) for some disordered mesoscopic cylinders when they 
are half-filled. The Hubbard interaction strength U is 
set at 1 and the impurity strength W is fixed to 2. In (a) 
the ground state energy level is shown for a cylinder with 
N = 5 (odd), while in (b) it is presented for a cylinder 
taking A^ = 8 (even). Quite interestingly we see that for 
the cylinder with odd N, the half flux-quantum period- 
icity of the lowest energy level disappears in the presence 
of impurity and it provides conventional 00 periodicity. 
Hence, for cylinders with odd N, 4>o/2 flux-quantum pe- 
riodicity will be observed only when they are free from 
any impurity. For the disordered cylinder with even N 
{N = 8), the lowest energy level as usual provides (f>o pe- 
riodicity similar to the impurity free cylinders containing 
even N. Apart from this periodic nature, impurities play 
another signiflcant role in the determination of the slope 
of the energy levels. The slope of the lowest energy level 
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FIG. 5: (Color online). Ground state energy level as a func- 
tion of flux (j) for half-filled disordered mesoscopic cylinders 
{M = 2) considering [/ = 1 and VF = 2. (a) iV = 5 and (b) 



decreases significantly compared to the perfect case, and 
therefore, a prominent change in current amplitude also 
takes place. 

To justify the above facts, in Fig.[S]we present the vari- 
ations of persistent currents with AB flux (j) for a half- 
flUed mesoscopic cylinder, described in the framework of 
NNH model, considering A^ = 15 and M ~ 2. The red 
curve represents the current for the ordered {W = 0) 
non-interacting {U = 0) cylinder. It shows saw-tooth like 
nature with flux (f> providing (j)o/2 flux-quantum period- 
icity. The situation becomes completely different when 
impurities are introduced in the cylinder as seen by the 
other two curves. The green curve represents the current 
for the case only when impurities are considered but the 
effect of electronic correlation is not taken into account. 
It shows a continuous like nature with ^o flux-quantum 
periodicity. The most important observation is that the 



current amplitude gets reduced enormously, even an or- 
der of magnitude, compared to the perfect cylinder. This 
is due to the localization of the energy eigenstates in the 
presence of impurity, which is the so-called Anderson lo- 
calization. Hence, a large difference exists in the current 
amplitudes of an ordered and disordered non-interacting 
cylinders and it was the main controversial issue among 
the theoretical and experimental predictions. Experi- 
mental verifications suggest that the measured current 
amplitude is quite comparable to the theoretical current 
amplitude obtained in a perfect system. To remove this 
controversy, as a first attempt, we include the effect of e- 
e correlation in the disordered cylinder described by the 
NNH model. The result is shown by the blue curve where 
U is fixed at 1.5. It is observed that the current am- 
plitude gets increased compared to the non-interacting 
disordered cylinder, though the increment is too small. 
Not only that the enhancement can take place only for 
small values of U, while for large enough U the current 
amplitude rather decreases. This phenomenon can be im- 
plemented as follows. For the non-interacting disordered 
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FIG. 6: (Color online). Persistent current as a function of flux 
(j> for a half-filled mesoscopic cylinder considering A'' = 15 and 
M — 2. The red line corresponds to the ordered case when 
[/ = 0, whereas the green and blue lines correspond to the 
disordered case {W = 2) when [/ = and 1.5, respectively. 

cylinder the probability of getting two opposite spin elec- 
trons becomes higher at the atomic sites where the site 
energies arc lower than the other sites since the electrons 
get pinned at the lower site energies to minimize the 
ground state energy, and this pinning of electrons be- 
comes increased with the rise of impurity strength W. 
As a result the mobility of electrons and hence the cur- 
rent amplitude gets reduced with the increase of impurity 
strength W. Now, if we introduce electronic correlation 
in the system then it tries to depin two opposite spin 
electrons those arc situated together due to the Coulomb 
repulsion. Therefore, the electronic mobility is enhanced 
which provides larger current amplitude. But, for large 
enough interaction strength, no electron can able to hop 
from one site to other at the half-filling since then each 
site is occupied either by an up or down spin electron 
which does not allow other electron of opposite spin due 
to the repulsive term U. Accordingly, the current ampli- 
tude gradually decreases with U. On the other hand, at 



less than half-filling though there is some finite probabil- 
ity to hop an electron from one site to the other available 
'empty' site but still it is very small. So, in brief, we 
can say that within the nearest-neighbor hopping (NNH) 
approximation electron-electron interaction does not pro- 
vide any significant contribution to enhance the current 
amplitude, and hence the controversy regarding the cur- 
rent amplitude still persists. 



C. Disordered cylinders with NNH and SNH 
integrals 

To overcome the existing situation regarding the cur- 
rent amplitude, in this sub-section, finally we make an at- 
tempt by incorporating the effect of second- neighbor hop- 
ping (SNH) integral in addition to the nearest-neighbor 
hopping (NNH) integral. 

A significant change in current amplitude takes place 
when we include the contribution of second- neighbor hop- 
ping (SNH) integral in addition to the NNH integral. As 
representative examples, in Fig. [7] we plot the current- 
flux characteristics for a half-filled mesoscopic cylinder 
considering TV = 15 and A/ = 2. The black, magenta 
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FIG. 7: (Color online). Persistent current as a function of 
flux (j) for a half-filled mesoscopic cylinder taking A'' = 15 and 
Af = 2 in the presence of NNH and SNH integrals. The black 
line corresponds to the ordered case when [/ = 0, whereas 
the magenta and gold lines correspond to the disordered case 
{W = 2) when f/ = and 1.5, respectively. Here SNH integral 
is fixed at —0.6. The currents shown by the red, green and 
blue lines for the ring described with NNH model (identical 
to Fig. [6| are re-plotted to judge the effect of SNH integral 
over NNH model much clearly. 

and gold lines correspond to the results in the presence 
of SNH integral, while the other three colored curves 
(red, green and blue) represent the currents in the ab- 
sence of SNH integral. Here we choose td = —0.6. The 
black curve refers to the persistent current for the perfect 
{W = 0) non-interacting ([/ = 0) cylinder and it achieves 
much higher amplitude compared to the NNH model (red 
curve). This additional contribution comes from the SNH 
integral since it allows electrons to hop further. In addi- 
tion it is also noticed that the current varies periodically 
with <j) providing 0q flux-quantum periodicity, instead of 




0o/2 as in the case of NNH integral model (red curve). 
Thus, it can be emphasized that (/)o/2 periodicity will be 
observed only when the cylinder is (a) free from impurity, 
(b) half-filled, (c) made with odd iV, and (d) described 
by the nearest-neighbor hopping model. The main focus 
of this sub-section is to interpret the combined effect of 
SNH integral and electron-electron correlation on the en- 
hancement of persistent current amplitude in disordered 
cylinder. To do this first we narrate the effect of SNH 
integral in disordered non-interacting cylinder. The na- 
ture of the current for this particular case is shown by the 
magenta curve of Fig. [71 It shows that the current am- 
plitude gets reduced compared to the perfect case (black 
line), which is expected, but the reduction of the current 
amplitude is quite small than the NNH integral model. 
This is due the fact that the SNH integral tries to de- 
localize the electronic states, and therefore, the mobility 
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FIG. 8: (Color online). Persistent current as a function of 
flux (j) for a half-filled mesoscopic cylinder taking A'' = 15 and 
Af = 2 in the presence of NNH and SNH integrals. The black 
line corresponds to the ordered case when [/ = 0, whereas 
the magenta and gold lines correspond to the disordered case 
{W = 2) when [/ = and 1.5, respectively. Here SNH integral 
is fixed at —0.8. The currents shown by the red, green and 
blue lines for the ring described with NNH model (identical 
to Fig. [6]) are re-plotted to judge the effect of SNH integral 
over NNH model much clearly. 

of the electrons is enriched. The situation becomes more 
interesting when we include the effect of Hubbard inter- 
action. The behavior of the current in the presence of 
interaction is plotted by the gold curve of Fig. [7] where 
we fix [/ = 1.5. Very interestingly we see that the current 
amplitude is enhanced significantly and quite comparable 
to that of the perfect cylinder. 

For better clarity of the results discussed above, in 
Fig. [8] we also present the similar feature of persistent 
current for other hopping strength of SNH integral. Here 
we set td = —0.8. From these curves we see that the 
current amplitude gets enhanced more as we increase the 
SNH strength. 

Thus, it can be predicted that the presence of SNH 
integral and Hubbard interaction can provide a persis- 
tent current which may be comparable to the measured 
current amplitudes. In this presentation we consider the 
effect of only SNH integral as a higher order hopping 
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integral in addition to the NNH model, and, illustrate 
how such a higher order hopping integral leads an im- 
portant role on the enhancement of current amplitude in 
presence of Hubbard correlation for disordered cylinders. 
Instead of considering only the SNH integral we can also 
take the contributions from all possible higher order hop- 
ping integrals with reduced hopping strengths. Since the 
strengths of other higher order hopping integrals arc too 
small, the contributions from these factors are reasonably 
small and they will not provide any significant change in 
the current amplitude. Finally, we can say that further 
studies are needed by incorporating all these factors. 



IV. CLOSING REMARKS 

To summarize, in the present work we have addressed 
the behavior of persistent current in an interacting nieso- 
scopic cylinder threaded by an Aharonov-Bohm flux 0. 
We have adopted a tight-binding Hamiltonian to describe 
the model quantum system and all the numerical calcu- 
lations have been done within a mean field approxima- 
tion. Using the generalized Hartree-Fock (HF) approx- 
imation, we have computed persistent current as func- 
tions of SNH integral, impurity strength W, AB flux (p, 
electron filling N^/N and system size N. Our numerical 
results have provided several interesting features and the 
present study may be helpful in understanding magnetic 
response in nano-scale loop geometries. 

The essential features observed from our analysis are 
as follows. 

(i) In the determination of the lowest energy level we see 
that the energy level varies periodically with cf) exhibiting 
both 00/2 and (f>o flux-quantum periodicities depending 



on the choices of the parameters describing the tight- 
binding Hamiltonian. 

(ii) In the NNH model, current amplitude gets signif- 
icantly reduced when impurities are introduced in the 
system. With the inclusion of Hubbard interaction (U), 
the current amplitude can be enhanced, though the en- 
hancement becomes too small compared to the experi- 
mental verifications. 

(iii) A significant change in the current amplitude takes 
place when the effect of SNH integral is taken into ac- 
count. The combined effect of SNH and Hubbard inter- 
action can provide the current which is quite comparable 
to the experimental realizations. So in short we can say 
that the conventional NNH model can be modified by 
incorporating the higher order hopping integrals. 

Throughout the analysis we have kept the width of the 
cylinders at a fixed value (M = 2), for the sake of our 
illustration. All these results are also valid for cylinders of 
larger widths. Here we have considered several important 
approximations by ignoring the effects of temperature, 
electron-phonon interaction, etc. Due to these factors, 
any scattering process that appears in the cylinder would 
have influence on electronic phases. At the end, we would 
like to say that we need further study in such systems by 
incorporating all these effects. 
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